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Abstract 

Using data collected during the 1992-1993 collider run at Fermilab, CDF has made 
measurements of correlated b quark cross sections where one b is detected from the 
lepton from semileptonic decay and the second b is detected with secondary vertex 
techniques. We report on measurements of the cross section as a function of the 
momentum of the second b and as a function of the azimuthal separation of the two 
b quarks, for transverse momentum of the initial b quark greater than 15 GeV. The 
vertex reconstruction techniques are valid over a large range in transverse momentum, 
starting at a minimum of 10 GeV. Results are compared to QCD predictions. 
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Studies of b production in pp collisions provide quantitative tests of perturbative QCD. 
For processes involving momentum transfers on the order of m^, the strong coupling constant, 
as, becomes relatively small and perturbative methods are expected to work well. Measure- 
ments of inclusive cross sections for pp — > bX have been made at CDF and UAl [Q. 
Consideration of the process pp —>■ bbX provides further opportunities for comparison of 
experiment and QCD calculations. 

We will make measurements of correlated b quark cross sections, where one b is detected 
from the lepton from semileptonic decay and the second b is detected with secondary vertex 
techniques. We report on measurements of the cross section as a function of the transverse 
momentum of the second b {dat/dET) and as a function of the azimuthal separation {da/d6(f)) 
of the two b quarks, for transverse momentum of the initial b quark greater than 15 GeV. 

This paper presents the first correlated cross sections measured at CDF using a combina- 
tion of lepton and vertex tags to identify b events. Starting with a well identified fi candidate, 
we look for the presence of a second b with secondary vertex techniques. Using jets with 
Et > 10 GeV, we measure the cross section as a function of the b transverse momentum and 
also the cross section as a function of the azimuthal angle between the b and b quarks. For 
purposes of notation, we will consider the fi as coming from a b quark, and the jet as coming 
from a b quark. 

The CDF has been described in detail elsewhere 0]. The tracking systems used for this 
analysis are the silicon vertex detector (SVX), the central tracking chamber (CTC), and the 
muon system. The SVX and CTC are located in a 1.4 T solenoidal magnetic field. The 
SVX consists of 4 layers of silicon-strip detectors with r — readout, including pulse height 
information with a total active length of 51 cm. The pitch between readout strips is 60 fim 
and a spatial resolution of 13 /im has been obtained. The first measurement plane is located 
2.9 cm from the interaction point, leading to an impact parameter resolution of ~ 15 /im for 
tracks with transverse momentum, p^, greater than 5 GeV/c. The CTC is a cylindrical drift 
chamber containing 84 layers, which are grouped into alternating axial and stereo superlayers 
containing 12 and 6 wires respectively, covering the radial range from 28 cm to 132 cm. The 
central muon system consists of two detector elements. The Central Muon chambers (CMU), 
located behind ~ 5 absortion lengths of material, provide muon identification over 85% of 
(f) for the pseudorapidity range \ri\ <0.6, where t] = — ln[tan(^/2)]. This rj region is further 
instrumented by the Central Muon Upgrade chambers (CMP), located after ^ 8 absorption 
lengths. The calorimeter systems used for this analysis are the central and plug systems. 
The central subtends the range \r]\ < 1.1 and spans 27r in azimuthal coverage. The plug 
subtends the range 1.1 < \t]\ < 2.4, again with 27r azimuthal coverage. 

CDF uses a three-level trigger system. At Level 1, muon candidate events are selected 
with a trigger that requires the presence of a hit pattern, consistent with p^ > 6 GeV/c, 
in the CMU chambers and confirming hits in the CMP chambers. At Level 2, the trigger 
requires that the CMU chamber track match a track found in the CTC, using the Central 
Fast Track processor, with p^ > 9.2 GeV/c. At Level 3, the trigger requires a muon track 
in both the CMU and CMP chambers matched with a track found in the CTC, using the 
offline track reconstruction algorithm, with p^ > 7.5 GeV/c. 

An inclusive muon sample is formed from the muon triggered sample with the following 
selection criteria: (1) reconstructed CTC track, pt > 9 GeV/c, (2) muon tracks exist in both 
the CMU chambers and CMP chambers, (3) the of the match betweent muon track and 
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the extrapolated CTC track match be less than 9 in the transverse view of both CMU and 
CMP chambers and be less than 12 for the longitudinal view in the CMU chambers, and 
(4) the track extrapolate and is found in the SVX fiducial region. There are 145784 events 
passing all requirements in this data sample. 

Monte Carlo samples for b and c quarks are produced using ISAJET version 6.43 The 
CLEO Monte Carlo program P| is used to model the decay of b hadrons. b quarks produced 
using the HERWIG Monte Carlo are also used for systematic studies. One of the final 
state b partons is required to have pt > 15 GeV. Events with a /i with p^ > 8 GeV are 
passed through the full CDF simulation and reconstruction package. The simulation used 
an average b lifetime of cr = 420 /im. 

The /i acceptance and efficiency has three parts to it: (1) the fiducial acceptance for 
muons coming from 6's with \r]\ < 1, (2) the fraction of 6's, > p™"-^ which decay to muons 
with pt > 9 GeV and (3) the trigger and identification efficiencies for 9 GeV muons. The 
first two factors have been studied using the Monte Carlo samples described above. The 
trigger and identification efficiencies for muons are studied from the data, using J/ip and Z° 
samples. All acceptance and efficiency numbers are summarized in table 0. 

The p^*"' value is chosen according to the standard interpretation, where 90% of muons 
with pf > 9 GeV come from b quarks with pt greater than p™". For this sample, the p]^^^ 
value is 15 GeV, from studying a sample of b quarks generated with pt > 10 GeV. 

We find that 18.6 ± 1% of muons with pt > 9 GeV pass through the CMU— CMP 
fiducial regions, where the error is statistical only. This value is for the entire r] distribution 
of produced b quarks. Restricting the |?7| of the parent b quark to be less than 1, we find 
39.2 ± 2% of muons with pj > 9 GeV pass through the CMU — CMP chambers. From a 
study of b quarks with pt > 15 GeV which decay to fi, we find that 10.7 ± 0.1 % (statistical 
error only) have pj > 9 GeV. We use the branching ratio B(6 ^ /i) = 0.108 ± 0.065 which 
comes from direct measurements at CLEO This fraction includes the sequential decay 
contribution, which has been scaled by the relative branching ratios of 6 /i to c ^ /i. 
We have varied the mean < z > used in the b fragmentation |p to study the effects in the 
geometric and kinematic acceptance. The systematic errors are therefore correlated. A la 
variation in the mean < z > gives a +1 — 6.5 % change in the geometric acceptance and a 
+9.7 — 10.8% change in the kinematic acceptance. 

The trigger efficiency is measured using independently triggered samples for each level 
of the system. The efficiency curves are then convoluted with the pt spectrum of muons, 
to get the efficiency for a muon with pt > 9 GeV. The combined efficiency of the LI, L2, 
and L3 triggers is measured to be 0.81 ± 0.024. The efficiency of the muon reconstruction 



algorithms has been studied in detail in the muon plus charm meson cross section work |10 
and is found to be 0.981 ± 0.003. Given the presence of reconstructed stubs in the chambers, 
the matching efficiency has been studied in a J/?/' sample and found to be 0.987 ± 0.013 for 
the matching cuts used in this analysis. The requirement that the CTC track be in the SVX 
fiducial region, in combination with the SVX tracking efficiency, has an efficiency of 0.68 ± 
0.021. 

The combined acceptance and efficiency for muons coming from b quarks with pt > 15 
GeV is 0.00239 +0.00030 - 0.00018. This correction will be applied to all cross section 
numbers presented below. Therefore, there is a common uncertainty of +12.4 — 7.5% on all 
the cross sections coming from the n acceptance and identification efficiency. 
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Geometric Acceptance 
Kinematic Acceptance 
Branching Ratio 


0.392 +0.005 -0.026 
0.107 +0.010 -0.012 
0.108 ± 0.065 


Trigger Efficiency 
Reconstruction Efficiency 
Matching Efficiency 
SVX requirement 


0.81 ± 0.024 
0.981 ± 0.003 
0.987 ± 0.013 
0.68 ± 0.021 


Combined Acceptance 
and Efficiency 


0.00239 +0.00030 - 0.00018 



Table 1: Summary of muon acceptance and efficiency numbers. The uncertainties in the 
kinematic and geometric acceptance are correlated. 



For the correlated cross sections, we look for additional jets in the sample, binning by 
jet Et and azimuthal opening angle. For each jet in the sample, we calculate a probability 
that the jet is consistent with coming from the primary vertex. Using input templates from 
b, c, and primary (light quark and gluon) jet Monte Carlo samples, we ffi the data as a 
sum of these three processes. Assuming that the presence of the /z candidate signals an 
independent b quark, we can then make correlated cross section measurements. Correcting 
for the acceptance and identification requirements and the b jet acceptance, and scaling 
by the integrated luminosity of the sample, we have differential cross section measurements. 

We make use of a jet probability algorithm |11| which asks the question "Is the ensemble 
of tracks in this jet consistent with being from the primary vertex?" This algorithm compares 
track impact parameters to measured resolution functions in order to calculate for each jet 
a probability that it is attached to the origin. This probability is uniformly distributed for 
light quark or gluon jets, but is very low for jets with displaced vertices from heavy flavor 
decay. The algorithm assigns a probability between and 1, where a probability value near 
1 means the jet is very consistent with being primary, and a probability value near means 
the jet is very inconsistent with being primary. 

For jet clusters identified in the calorimeter, the algorithm selects a set of tracks, p^ > 
1 GeV/c, within a cone of 0.4 around the jet axis to be used in the calculation of the jet 
probability. There are loose cuts to affiliate the tracks with the primary vertex, in addition 



to track quality requirements |Tl]]. We require that there be > 2 tracks passing the quality 
requirements for the calculation of the probability. 

Jets in this sample are required to have Ey > 10 GeV in a cone of radius 0.4, \ri\ < 
1.5, have at least 2 good tracks, and be separated from the muon in r/ — space by AR > 
1.0. The cone of radius 1.0 was chosen so that the tracks clustered around the jet axis were 
separated from the fi direction. All jet energies in this paper are measured energies, not 
including corrections for known detector effects jT^ . 

The b jet acceptance combines the fiducial acceptance of the SVX and the CTC, the 
track reconstruction efficiency, and fragmentation effects. The simulation is used to calculate 
the combination of these pieces. The acceptance represents the fraction of b quarks which 
produce jets with Et > 10 GeV, \ri\ < 1.5 and at least 2 good tracks inside a cone of 0.4 
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Range 


Acceptance 


10 


— 15 


32.9 ± 1.9 % 


15 


— 20 


46.3 ± 2.2 % 


20 


— 25 


48.8 ± 2.9 % 


25 


— 30 


50.8 ± 3.9 % 


30 


— 40 


49.7 ± 4.2 % 


40 


— 50 


49.8 ± 7.3 % 



Table 2: b jet acceptance as a function of jet for \ri\ < 1.5 (statistical errors only). 



S(f) Range 


Acceptance 


— 


22.5° 


7.25 ± 2.2 % 


22.5° 


— 45° 


4.52 ± 2.0 % 


45° - 


-67.5° 


36.9 ± 3.0 % 


67.5° 


— 90° 


51.4 ± 0.8 % 


90° — 


112.5° 


51.4 ± 0.8 % 


112.5° 


— 135° 


51.4 ± 0.8 % 


135° - 


- 157.5° 


51.4 ± 0.8 % 


157.5° 


— 180° 


51.4 ± 0.8 % 



Table 3: 5 jet acceptance as a function of 6(j) of the two b quarks, for Ej' > 10 and |?7| < 1.5 
(statistical errors only). The first three bins have lower acceptance due to the AR cut for 
identification of the two b jets. 



around the jet axis, where there is also a b quark which decays to a with pt > 9 GeV within 
the CMU-CMP acceptance. The b jet acceptance is calculated separately as a function of 
the jet Et and azimuthal opening angle between the two quarks. 

The average acceptance for the 6 is ~ 40%. It ranges from 32.9 ± 1.9% (statistical error 
only) for 10 < Et < 15 GeV to 49.8 ± 7.3% for 40 < Et < 50 GeV. For 50 < 22.5°, 
the acceptance is 7.3 ± 2.2%, while for 157.5° < 50 < 180°, the acceptance is 51.4 ± 0.8%. 
Tables ^ and ^ show the bin by bin values used in the differential cross section measurements. 

We have compared the values for the b jet acceptance from ISAJET samples to the 
acceptance from HERWIG samples. The acceptance agrees within the statistical error in 
the samples as a function of E^, differing at the 5% level. We will take this as an additional 
systematic uncertainty on the acceptance. In combination with a 10% uncertainty due to 
the vertex distribution for events in the SVX fiducial volume, we have a common 11.2% 
systematic uncertainty in all the jet acceptance numbers. 

We use a binned maximum likelihood fit |T3] to distinguish the b, c, and primary jet 



contributions in the sample. In this case, we use the results of the jet probability algorithm 
to define the variable that is used in the fit. We find that the logio(jet probability) shows 
stronger differentiation between b, c, and primary jets (see figure IT]) than the jet probability 
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and will use this variable in the fitting algorithm. We fit over the range -10 — in logio(jet 
probability), where the b, c, and primary contributions are constrained to be positive, but 
no other constraints are included in the fit. We use smoothed Monte Carlo distributions 
from b and c samples as the input shapes for heavy flavor to decrease the effects of limited 
Monte Carlo statistics. We model the primary jets with an exponential distribution, since a 
logarithm transforms a uniform distribution to an exponential distribution. 

We have explored the effect of different Monte Carlo samples to form the input shape used 
in the fit. Using different input b Monte Carlo samples compared to a test distribution made 
with independent Monte Carlo samples shows a 5% change in the fit fractions. Changing the 



average b lifetime by 6% [|T4[ changed the fit fraction by 3%. We include a 5.8% systematic 



uncertainty to our fit results to account for systematic errors in the fitting procedure and 
uncertainty of the b lifetime. 

In figure |^, we show the distribution of logio(jet probability) for all jets in the /i sample, 
overlayed with the fit results. In this sample, the fit predicts 2620 ±97 6 jets, 2085 ± 180 c 
jets, and 13103 ± 161 primary jets for a total of 17808. There are 17810 events in the data 
sample. Figure ^ shows three comparisons of the data and fit results, showing the bin-by-bin 
difference in the results, the bin-by-bin difference divided by the errors, and the distribution 
of the difference divided by the errors. In these distributions, the errors are a combination 
of the statistical errors in the data points and the overall scale error in the fitted values. We 
do not include any error on the Monte Carlo shapes. From these distributions, we can see 
that the inputs model the data well. The difference divided by the errors has a mean of -0.04 
and RMS of 1.08. 

We have applied the fitting algorithm to several different samples of jets, where we have 
made different requirements on the transverse energy of the jet or the azimuthal separation 
between the /i and the jet. Applying the acceptance and efficiency numbers for the /i can- 
didate, the b jet acceptance and dividing by the integrated luminosity of the sample, we 
convert the number of fit b jets into differential cross sections. All of the cross section num- 
bers presented have a set of uncertainties in common. These come from the n acceptance and 
identification (+12.4 — 7.5%), the fitting procedure (5.8%), the vertex acceptance (10%), 
the b jet acceptance (5%), and the luminosity normalization (3.6%). The total common 
systematic error is +18.3 — 15.4%. 

We look at the distribution of the jet in the event, using 6 bins to cover the region 
between 10 GeV and 50 GeV in transverse energy. In each bin, we do an independent 
fit of the logio(jet probability) distribution and then correct for the acceptance. Table § 
contains a summary of the number of b jets, number of total jets, and the cross section in 
each Et bin considered. Figure ^ shows the values of the cross section as a function of jet 
Et for yu with > 9 GeV. 

We also measure the cross section as a function of the azimuthal separation between the 
jet and the /x. The measured jet direction is a good measure of the initial b direction, while 
the n follows the b direction. We have chosen to use bins of 22.5° in 6(j), which is the width 
of the spread between the b direction and the n direction. With this bin size, the azimuthal 
separation of the jet and yU is a good approximation of the azimuthal separation of the two 
quarks. Table ^ contains a summary of the number of jets, number of fitted b jets, and the 
cross section in each 50 bin considered. 

Figure |^ shows a comparison to a theoretical prediction, using the Mangano-Nason-Ridolfi 



6 



Et range 


Number of Jets 


Number of fit 6 jets 


Cross Section (nb/GeV) 


10 — 15 


5695 


650 ± 52 


5.49 ± 0.63 


15 — 20 


4083 


647 ± 48 


3.88 ± 0.41 


20 — 25 


2680 


452 ± 40 


2.57 ± 0.32 


25 — 30 


1774 


323 ± 33 


1.77 ± 0.25 


30 — 40 


1915 


308 ± 32 


0.86 ± 0.13 


40 — 50 


835 


145 ± 21 


0.41 ± 0.085 



Table 4: Cross section summary as a function of jet E^. There is a common systematic 
uncertainty of +18.3 — 15.4% not included in the experimental points. 



50 Range 


Number of Jets 


Number of fit b jets 


Cross Section (nb/22.5°) 


— 22.5° 


83 


11.7 ± 7.4 


2.24 ± 1.58 


22.5° — 45 ° 


143 


35.8 ± 10.6 


10.88 ± 5.78 


45° — 67.5° 


364 


73 ± 16 


2.25 ± 0.66 


67.5° — 90° 


442 


81 ± 17 


2.19 ± 0.48 


90° — 112.5° 


624 


122 ± 20 


3.30 ± 0.58 


112.5° — 135° 


1188 


195 ± 27 


5.27 ± 0.80 


135° — 157.5° 


3394 


477 ± 43 


12.9 ± 1.40 


157.5° — 180° 


11752 


1624 ± 76 


43.9 ± 3.34 



Table 5: Cross section summary as a function of 6(j) between the jet and jj,. There is a 
common systematic uncertainty of +18.3 — 15.4% not included in the experimental points. 
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calculation |]T5| with nib = 4.75 GeV, MRSB structure functions, and /i^ = ml + < pt for 
\ri^\ < 1-, Pf > 15 GeV, \ri^\ < 1.5, and > 10 GeV as input. We have smeared the parton pt 



distribution with the expected energy scale and resolution of the CDF detector [|T^ , and have 



normalized the 50 distribution to the cross section with smeared pt > 10 GeV. The shapes 
of the theoretical prediction and the experimental data agree well, especially for S(f) > 90°. 
Note that there is a large change in the acceptance for 50 < 60° (see table H) due to the AR 
separation requirement on the n jet system. 

In conclusion, we have presented the first measurements of correlated bb cross sections 
using a combination of lepton and vertexing techniques to identify the two 6's. We mea- 
sure the differential transverse energy cross section for the b, dcr^/dE^, and the differential 
azimuthal cross section between the two quarks, da/d6(f). The shape in the 6(f) distribution 
agrees well with theoretical predictions, but the overall normalization is roughly a factor of 
1.3 higher than predicted. 
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National Science Foundation; the Italian Istituto Nazionale di Fisica Nucleare; the Ministry 
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Figure 1: The logioQet probability) distributions used as inputs to the fitting program. The 
b and c shapes are smoothed versions of Monte Carlo distributions, while the primary shape 
is an exponential function. The three distributions are normalized to equal area and shown 
on the same vertical scale. 
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Figure 2: For 10 GeV jets in the // sample, we show the data distribution overlayed with the 
fit results. Statistical errors on the data and scale errors on the fit results are included. The 
fit results model the data well over the entire range of the fit. 
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Figure 3: Various comparisons of the data distribution and the fit results. We show the bin 
by bin difference between the data and the fit resuhs, the bin-by-bin difference scaled to the 
errors, and the distribution of the difference scaled to the errors. In all cases, the errors are 
the statistical error in the data points plus the overall scale error in the fitted values. 
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Figure 4: The distribution of the second 6 jet, > 10 GeV and \r]\ < 1.5, cross section as 
a function of jet E^, given a /i with Pt > 9 GeV present in the event. There is a common 
systematic uncertainty of +18.3 — 15.4% not included in the experimental points. 
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Figure 5: The distribution of the cross section as a function of the azimuthal angle between 
the fi and the jet, with Ej- > 10 GeV and |?7| < 1.5, overlayed with a theoretical prediction. 
The theoretical prediction has been normalized to the expected cross section after including 
the effects of detector smearing. There is a common systematic uncertainty of +18.3 — 15.4% 
not included in the experimental points. Note that the acceptance changes rapidly and has 
large uncertainty for the region 50 < 60° due to the effects of the AR separation cut. 
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